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We theoretically analyze a discrete Schrodinger chain with hopping to the first and second neigh¬ 
bors, as can be realized with zigzag arrangements of optical waveguides or lattice sites for cold 
atoms. Already at moderate values, the second-neighbor hopping has a strong impact on the band 
structure, leading to the emergence of a new extremum located inside the band, accompanied by 
a van Hove singularity in the density of states. The energy band is then divided into a subcritical 
regime with the usual unique correspondence between wave number and energy of the travelling 
waves, and a supercritical regime, in which waves of different wave number are degenerate in energy. 
We study the consequences of these features in a scattering setup, introducing a defect that locally 
breaks the translational invariance. The notion of a local probability current is generalized beyond 
the nearest-neighbor approximation and bound states with energies outside the band are discussed. 
At subcritical energies inside the band, an evanescent mode coexists with the travelling plane wave, 
giving rise to resonance phenomena in scattering. At weak coupling to the defect, we identify a 
prototypical Fano-Feshbach resonance of tunable shape and provide analytical expressions for its 
profile parameters. At supercritical energies, we observe coupling of the degenerate travelling waves, 
leading to an intricate wave packet fragmentation dynamics. The corresponding branching ratios 
are analyzed. 

PACS numbers: 03.75.Lm, 67.85.Hj, 42.79.Gn 


I. INTRODUCTION 

In many fields of physics, understanding and control¬ 
ling the propagation and scattering properties of waves 
is of profound importance [T]. Matter wave scatter¬ 
ing from a localized defect within the framework of the 
Schrodinger equation is one of the most basic, yet at the 
same time fundamentally important, problems of single 
particle quantum mechanics. This is particularly true in 
view of present-day technology which makes it possible to 
fabricate devices that are adequately described by basic 
quantum mechanical models, for example in electronics 
El [3]. At the same time, ultracold atom experiments 
are also entering the regime of mesoscopic physics mi, 
promising to make the enormous toolbox for cold atomic 
gases accessible in quantum transport investigations. 
The complexity of a scattering problem is not solely de¬ 
termined by the actual scatterer, but crucially depends 
on the structure of the asymptotic regions supporting the 
in- and outgoing waves. For instance, the interplay of dif¬ 
ferent transverse modes of quasi-one dimensional quan¬ 
tum waveguides gives rise to a rich variety of resonance 
phenomena, even when considering only a point-like scat¬ 
terer HHII], see also the confinement induced resonances 
in ultracold two-body collisions [laiis]- Nontrivial ef¬ 
fects can also arise from an intricate discrete structure 
of the asymptotic leads, e.g. in carbon nanotubes or 
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graphene [HHisi- Here, the detailed structure of the 
leads can often be modeled theoretically within a tight- 
binding approximation m, representing the material by 
an abstract multiply-connected lattice of discrete sites. 
At the same time, tight-binding lattices beyond the stan¬ 
dard first-neighbor approximation have also received con¬ 
siderable recent attention in the realm of cold atoms. 
Protocols employing ions or Rydberg atoms have been 
proposed for the simulation of bosonic Hubbard mod¬ 
els with long-range hopping [iHHin]- A powerful method 
for enhancing specifically the second-neighbor hopping 
term in a Hubbard chain relies on rearranging the lat¬ 
tice sites in a zigzag geometry m- This has been imple¬ 
mented E2] in evanescently coupled optical waveguide ar¬ 
rays ES], whose versatile applicability for the simulation 
of quantum scattering phenomena was also demonstrated 
in the realization of Klein tunneling EH ES] • Effects 
of the second-neighbor hopping on nonlinear excitations 
EIlEiET], Bloch oscillations EHHSO], wave localization 
m, or the Mott insulator transition E2] have been in¬ 
vestigated. Recently, several studies have explicitly ad¬ 
dressed ultracold atoms in zigzag lattice geometries, both 
with interactions EH El] and without ESI- Here, the 
zigzag geometry could be obtained experimentally by sin¬ 
gling out a strip from an extended two-dimensional trian¬ 
gular optical lattice ESj? oi* using recently demonstrated 
techniques for designing essentially arbitrary in-plane op¬ 
tical potentials EH EE]. 

The hallmark signature of sizable second-neighbor hop¬ 
ping in a noninteracting tight-binding chain is a defor¬ 
mation of the band structure, causing the emergence of 
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a split band edge m (see also m) and a correspond¬ 
ing van Hove singularity in the density of states located 
inside the band [50]. We will investigate here the in¬ 
fluence these features have when the zigzag lattice con¬ 
stitutes the asymptotic region of a scattering problem. 
This setup has the advantage of providing nontrivial, 
yet generic, band structure features, while at the same 
time being accessible to analytical methods. In fact, 
for similar reasons, it has been briefly suggested before 
as a minimal model for illustrating the Fano-Feshbach 
resonances gulls in elastic light scattering off an ob¬ 
stacle m, but only a restricted parameter range of the 
model was considered there. Beyond-flrst neighbor dis¬ 
crete Schrodinger lattice models also arise in the descrip¬ 
tion of macromolecules such as DNA [33], and in this con¬ 
text the zigzag defect problem has been discussed in m, 
with a focus, however, on bound states instead of scatter¬ 
ing properties. Going beyond these works, we provide a 
comprehensive analysis of scattering from a localized de¬ 
fect in the zigzag geometry for a wide range of lattice pa¬ 
rameters and energies, identifying pronounced resonance 
features and the possibility of wave packet fragmentation 
induced by the second-neighbor hopping. 

Our presentation is structured as follows: In Sec. we 
introduce the model and give a detailed discussion of the 
band structure features in the absence of the defect. Sec. 
|III| explores the constraints imposed on stationary scat¬ 
tering solutions by continuity. In Sec. El we discuss 
bound states at the defect, before entering the discus¬ 
sion of scattering states. Here, we distinguish two ener¬ 
getic regimes, a subcritical one in which the interplay of 
a closed and an open channel induces a Fano-Feshbach 
resonance (Sec. |v|), and a supercritical one in which two 
open channels coex ist a nd can be coupled to each other 
via scattering (Sec. VI). In both regimes we provide ana¬ 
lytical results for the stationary scattering solutions and 
show corresponding simulations of the wave packet dy¬ 
namics. Finally, we summarize and conclude in Sec. m 
and give a brief outlook on future perspectives. 


II. SETUP AND BAND STRUCTURE 

The discrete Schrodinger model we will be concerned 
with in the following has the general form 

idripj =- {ttjlpj+l + 

- {4,ji’j+2 + t2,j'>Pj-2) + VV’i • (1) 

Here, r denotes time, j G Z is the site index, and we al¬ 
low for site-dependent first {t^j) and second {t^j) neigh¬ 
bor hopping and an on-site potential term Vj. The hop¬ 
ping matrix is constrained by '^t,j — ^ 2 ,j+ 2 ' 

We use dimensionless units throughout. Such discrete 
Schrodinger models arise in a variety of contexts m, 
with correspondingly different physical interpretations of 
the terms. When thinking of ultracold atoms, Eq. 0 
describes the dynamics of a single particle (or a conden¬ 
sate of noninteracting bosons m) in a lattice poten¬ 


tial within the lowest-band approximation, i.e. taking 
into account a single localized Wannier mode per lat¬ 
tice site. We will not consider interactions here, not¬ 
ing that for many atomic species an effectively nonin¬ 
teracting limit can be prepared experimentally making 
use of internal-state Feshbach resonances m- Moreover, 
for the nearest-neighbor chain scattering results from 
the noninteracting system have been found to provide 
valuable information even in the presence of interactions 

[3HI SSI- 

Let us first discuss the homogeneous system with Vj = 
0, tfj = ti >0, = ^2 > 0 for all j. A sketch of this 

setup is shown in Fig. illustrating also the connection 
to the zigzag arrangement of lattice sites that this model 
is tailored for. In the zigzag lattice, the opening angle 0 
effectively tunes the relative distances between first and 
second (in index) neighbor sites, and thereby also the rel¬ 
ative values of the hopping amplitudes to the first and 
second neighbors ED ESI- Hopping matrix elements be¬ 
yond the second neighbor are neglected. 



Figure 1. (Color online) Sketch of the homogeneous discrete 
Schrodinger chain with hopping to the first and second neigh¬ 
bor and the zigzag arrangement of lattice sites which provides 
a realization of this model. Here, the opening angle 0 tunes 
the relative strength of t 2 compared to ti . 

In the homogeneous case, the model is invariant under 
spatial translations and Eq. 0 admits stationary plane 
wave solutions of the form ipj oc exp [i(/cj — E(/c)r)], 
where k denotes the quasi-momentum which can be re¬ 
stricted to the first Brillouin zone —tt < k < it. The 
corresponding dispersion relation is given by 

E{k) = —2ti cos k — 2^2 cos 2k. ( 2 ) 

This dispersion curve is shown in Eig. [^a) for the case 
ti = ^2 = 1- Characteristically, when increasing the 
second-neighbor hopping t 2 the global minimum of E{k) 
remains at /c = 0, Eq = E{0) = —2{ti + ^ 2 ), while the 
maximum of the dispersion is shifted away from k = it 
(which itself transforms into a local minimum) as soon 
as t 2 > ti/4, which is the case we will mostly focus on 
in the following. In particular, this implies that above 
a certain energy Ec within the band the dispersion is 
degenerate, in the sense that for E > Ec there are two 
distinct quasi-momenta ki^k 2 with \ki\ 7 ^ \k 2 \ satisfying 
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E{ki) = E{k 2 ). Remarkably, the group velocity 
diE 

^ + 4^2 sin 2k (3) 

dk 

is thus no longer in unique correspondence to the energy: 
The zigzag model permits wave packets of the same en¬ 
ergy that travel at different velocities in the lattice. The 
group velocity curve v{k) is also shown in Fig. [^a). 



Figure 2. (Color online) (a) Dispersion relation E{k) (solid 
red line) and group velocity v{k) (dashed grey line, shading 
highlights the intervals of positive v) and (b) normalized den¬ 
sity of states for the homogeneous zigzag lattice model with 
ti = t2 = 1. 


waves are coupled to each other. The most fundamen¬ 
tal framework to study this is a scattering setup, with a 
localized defect that breaks the homogeneity, while the 
asymptotic zigzag regions are undisturbed. Specifically, 
we will consider a localized on-site potential, Vj = Vdj^Q. 
Away from this defect, in the asymptotic regions of the 
lattice, we again assume first- and second-neighbor hop¬ 
pings > 0 and ^2 > 0 which are independent of the 
site index. In general, a local manipulation of the lattice 
that causes the on-site potential shift at the defect site 
will also affect the hopping matrix elements to this site, 
cf. the discussion in m- To keep the number of pa¬ 
rameters tractable, we will partially account for this by 
assuming that first- and second-neighbor hoppings that 
connect to the defect site are rescaled compared to the 
respective background values by a common factor 7 > 0. 
A local tuning of the hopping parameters only (without 
significantly altering the on-site potentials) is possible 
via a local variation of the inter-site distances [5T]. We 
will thus assume 7 and V to be essentially independently 
tunable parameters. The scattering setup is sketched in 
Fig.|3l 


For our model, the critical energy above which this 
travelling wave degeneracy occurs (given ^2 > ti/A) ex¬ 
plicitly reads as 


Ec = E{±kc) = 2{ti - 12 ), kc = acos (^1 - , 

while the split upper band edge is located at 

Efn — E{Ekjyi) = 2^2 T ^^5 ~ acos ( 4^ ) * 


In the limit of dominant second-neighbor hopping, ^2 ^ 
ti, this implies /Cc ^ 0, ^ 7r/2, such that for large ^2 

eventually the degeneracy region covers the whole band. 
In contrast, for weak second-neighbor hopping, ^2 < ti/4, 
Ec = 2(ti — ^2) marks the upper band edge dX k = Ttt, 
while Em = 2^2 + t^/(4t2) lies outside the band. 
Consider ^2 > ti/4, then the qualitative deformation of 
the equienergy “surface” in /c-space when crossing the 
critical energy Ec may be thought of as the basic ingre¬ 
dient for observing a topological transition in the cor¬ 
responding fermionic many-body system in the sense of 
m- For a Fermi energy near E^ the Fermi surface topol¬ 
ogy will be sensitive to small variations of the model pa¬ 
rameters. Characteristically, the formation of additional 
stationary points of the dispersion curve at the critical 
energy Ec leads to the emergence of a van Hove singular¬ 
ity in the density of states which is shown in Fig. |^b). 
The above considerations already suggest that by tuning 
the hopping parameters one can control the propagation 
and dispersion properties of wave packets in the homoge¬ 
neous zigzag lattice, as has been suggested in m- Here, 
we analyze the effects of the band structure deformation 
induced by the second-neighbor hopping if the transla¬ 
tional invariance of the system is broken and the plane 



Figure 3. (Color online) Sketch of the zigzag lattice containing 
a defect site with on-site potential V. Hopping to this defect 
site is modulated by a factor 7 , compared to the asymptotic 
homogeneous hoppings. 


Then, after the separation 'ipjir) = (j)j exp(—iF^r) with 
a time-independent (j)j and the energy E, the stationary 
discrete Schrodinger system of equations explicitly reads 

E(pj = —ti {(pj-i + — h (0J-2 + 0 ^+ 2 ) (4) 

for \j\ > 2, while 

E<^±2 = —ti (0±i + ^±3) ~ h (7^0 + ^±4) 5 
E(f)±l = —ti (700 + 0±2) — h (0IF1 + 0±3) , 

{E - H)0o = -7^1 (0-1 + 0l) - 7^2 (0-2 + 02) . 

(5) 


Let us first discuss the asymptotic regions of \j\ > 2, gov¬ 
erned by Eq. Here, the homogeneous zigzag lattice 
discussed above is recovered. However, since we now look 
for piecewise solutions on the semiaxes j > 2, j < —2 
only, we cannot restrict to the travelling plane waves, 
but need to take into account the possibility of evanes¬ 
cent waves which exponentially decay for j ±00, re¬ 
spectively. The fundamental solutions in the asymptotic 
regions are thus of the form 0^ oc ex.p{iKj) with a com¬ 
plex K. From Eq. K needs to satisfy the dispersion 
relation E = E{K) as in Eq. (§. Eor our model, we can 
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explicitly invert this to K{E)^ with the result 


cos Ki^ 2 {E) 




E 1 


(6) 


where for definiteness we choose the upper sign for Ki, 
the lower for K 2 . At a given energy E the equation 
E{K) = E thus has four complex solutions coming in 
pairs ±i^2- These solutions, again for ti = ^2 = 1, 
are visualized in Fig. at different values of E. Qualita¬ 
tively, the results do not change for other values of the 
hopping parameters as long as t 2 > ti/4. We can dis¬ 
tinguish four different regimes. For E < Eq^ i.e. below 
the lower band edge, we find one pair of solutions that 
is purely imaginary, while the other has a non-vanishing 
imaginary part and a constant real part of tt. Corre¬ 
spondingly, in this interval we have found no travelling 
waves. At Eq, the former pair of solutions turns purely 
real: There is now one pair of real solutions (correspond¬ 
ing to the travelling plane waves found before) and one 
pair of evanescent solutions of the form EK 2 = tt + m 
with K G M.. At Ec, the imaginary part of this second 
pair also goes to zero, and for Ec < E < E^ we are 
left with two real solutions, corresponding to the two de¬ 
generate travelling waves we have observed in the band 
structure. Finally, at Em there is a pairwise collision of 
the solutions in the complex plane upon which they leave 
the real axis and form a complex quartet. 
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Figure 4. (Color online) Complex solutions K(E) according 
to Eq. (|^ for ti = t 2 = 1 at different values of the energy 
E. Values oi Ki are indicated as red circles, values of K 2 
by grey squares, while arrows indicate the direction in which 
these solutions move with increasing E. Qualitatively, the 
same picture is found for any t 2 > ti/4. 


Let us contrast this to the case of weak second- 


neighbor hopping, t 2 < ti/4. Then below Ec (which in 
this case denotes the upper band edge) the situation is 
qualitatively the same as in panels (a) and (b) of Fig. 
with the Ki pair of solutions turning real at Eq and then 
moving towards Ttt with increasing energy. However, 
now at Ec it is not the K 2 pair of solutions that also 
reaches the real axis, but instead the Ki solutions reach 
±7r first and branch off along the imaginary axis, forming 
a second pair of the form ±(7r + in) for Ec < E < Em- 
Eventually, at Em there is again a pairwise collision (now 
taking place at the axes of Ttt real part) upon which a 
complex quartet in the plane is formed. For ^2^0 
this second collision is pushed to infinite energies and 
the nearest-neighbor scenario is recovered. 

From these considerations, one can already estimate 
the gross overall scattering features in the different 
energy ranges for ^2 > ti/4: Outside the band, below 
Eq or above Em^ there are no travelling wave solutions 
at all, but bound states at the defect are expected. In 
the subcritial energy regime, i.e. for Eq < E < E^ 
there is one travelling wave (open scattering channel) 
coexisting with an evanescent mode (closed channel), 
which may give rise to scattering resonances. Finally, for 
supercritical energies with Ec < E < Em^ both channels 
are open and may be mixed during scattering events. 
These three different regimes will be analyzed in detail 
in the following. 


III. CONTINUITY EQUATION 

Before entering into the discussion of the scattering 
properties of the point defect, let us first discuss the con¬ 
straints on the scattering coefficients due to local prob¬ 
ability conservation. Globally, the discrete Schrodinger 
equation Q preserves the norm A local conti¬ 

nuity equation is obtained by multiplying Eq. Q by 7/;* 
and taking the imaginary part, which results in 

(7) 

This identity holds at each j and independently of the 
set of on-site potentials {Vj}. It may be thought of as 
a Kirchhoff-type balance equation, relating the change 
of probability at a site to the currents in the four links 
connected to this site. Now, in spite of the multi¬ 
connectedness of the lattice due to the second-neighbor 
hopping, we can rewrite this in the form of a standard 
one-dimensional continuity equation by defining a local 
current as 


2 a-l 

■ (8) 

a=lp=0 

Then Eq. 0 reads as 

dr\lpj\‘^ + Jj — Jj-1 = 0 , 


(9) 
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where the symmetry j = t^x.j+a hopping ma¬ 

trix has been used. In particular, Eq. implies that 
for a stationary solution the current is a global constant, 
Jj = J for all j. In scattering scenarios, this equation re¬ 
lates the coefficients of the fundamental solutions in the 
asymptotic regions left and right of the defect [ 52 ], see 
below. We note that a generalization of the above con¬ 
cepts, in particular of the expression for the local current, 
for one-dimensional lattices with more extended hopping 
(beyond the second neighbor) is possible. 


IV. BOUND STATES 


We now turn to the zigzag scattering setup as shown in 
Fig-i In this section we briefly discuss bound states at 
the defect, having energies E outside the band. Then, as 
seen above, all solutions of E = E{K) have a nonvanish¬ 
ing imaginary part. Let Ai, A2 denote the two solutions 
whose imaginary part is positive (to prevent asymptotic 
exponential growth), then we search for bound states in 
the form 


-iK2j 


pj = Aie * -h A2e 

pj = 


J < 0 , 
j > 0, 


with complex coefficients Ai,A2,Bi,B2. Inserting this 
into Eqs. (§ yields a homogeneous system of linear equa¬ 
tions for (Ai, A2, 5 i, ^2, 0o) which has nonzero solutions 
only if the coefficients satisfy 


Bi — Ai, B2 — A2 
sin Ki 


sin Ki 
sin A2 


^ 1 , 


700 = 1 - 


sin Ao 


Ai. 


( 10 ) 

( 11 ) 


and simultaneously 

^ = £’(1-7^) 

2n»sing, sin 

sm Ko — sm iv 1 V ^ 


+ - 


( 12 ) 


Both below and above the band it is readily found that 
the expression on the right hand side of Eq. ( 12 ) is real. 


Thus, one can conclude that for any bound state energy 
El) outside the band and any fixed 7, there is precisely 
one value of the defect potential V that produces pre¬ 
cisely one bound state at this energy. This V{Ei),j) is 
given by Eq. Fig.ii: a) shows the corresponding 

dependence of V (horizontal axis) on the bound state 
energy E5 (vertical axis) for different values of 7. Con¬ 
versely, it illustrates that for 7 = 1 each defect potential 
V 0 also produces precisely one bound state (which 
lies below the band for V < 0 and above the band for 
U > 0 ). In contrast, if 7 < 1 then for U in a re¬ 
gion around zero there is no bound state at all, while 
if 7 > 1, for V in a region around zero there are two 
bound states (one above, one below the band). Eor large 
values of the defect potential |U|, the bound state en¬ 
ergy El) ^ V, irrespectively of 7, as is expected from 


perturbation theory. This limit is accompanied by an 
increasing localization of the eigenmode when energeti¬ 
cally moving away from the band. Obviously, also in the 
decoupling limit 7 ^ 0 we find Eij ^ V and asymp¬ 
totically perfect localization at t he defe ct. The relations 
between the coefficients in Eqs. dioln] ) ensure that the 
global phase of the wave function can be chosen such that 
it is real (as is to be expected from time-reversal symme¬ 
try). Thus, the continuity equation is trivially satis¬ 
fied with J = 0. Qualitatively, these features have been 
observed before for models with nearest-neighbor hop¬ 
ping only m and for weak second-neighbor hopping m- 
In the limit of ^2 ^ 0 , the imaginary part of K2 diverges, 
such that sin Ai/ sin A2 ^ 0 . Then Eq. reduces to 
V = — 2 ti cos + 2ti7^ exp(iAi), from which it follows 
that exp(—iAi) = —V// (2ti)^ + 27^ — 1, re¬ 
covering the nearest-neighbour result of [ 5 T], while at 
the same time Eqs. imply that ^2,^2 ^ 0, 

700 ^ Ai. 




Figure 5. (Color online) (a) Energy of bound state(s) vs. 
defect potential for varying coupling strength 7 and ti — t 2 — 
1. The energy band in which propagating solutions exist is 
indicated in grey. (b,c) Profiles of the two bound states (above 
and below the band) atU = ± 2,7 = 1.5, as marked by brown 
circles and dashed vertical lines in (a). 


A crucial difference to the nearest-neighbor model, 
however, is the absence of symmetry between bound 
states above and below the band. When ^2 = 0 , the 
staggering transformation (pj {—lycpj maps a station¬ 
ary solution at energy E for the potential strength V to 
a solution at energy —E for the potential strength —V. 
In particular, this implies that the bound state(s) at EV 
and —V share the same density profiles and differ only 
in their local phases, cf. m- Invariance under the stag¬ 
gering transformation does not hold in the presence of 
the second-neighbor hopping term. Thus, upon chang¬ 
ing the sign of the defect potential U, the corresponding 
bound states feature different density profiles. This can 
be seen in panels (b) and (c) of Eig. Notably, inter- 
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ference between the two evanescent waves constituting 
the bound state allows for a nonmonotonic density decay 
away from the defect site. This feature is most apparent 
in the extreme limit of t 2 ^ ti, where the zigzag lattice 
separates into two essentially decoupled chains and thus 
the bound states at the defect will dominantly populate 
every second site. 


V. SCATTERING AT SUBCRITICAL 
ENERGIES 


From here on we focus on energies lying inside the 
band of propagating states. In this section we discuss 
scattering off the defect at subcritical energies satis¬ 
fying Eq < E < Ec. For weak second-neighbor hopping, 
h < ^i/4, this energetic regime covers the whole energy 
band {Ec coincides with the upper band edge), while for 
h > ^i/4 the critical Ec lies inside the band. The follow¬ 
ing discussion applies to both cases likewise. 

For Eq < E < Ec^ there is a unique real wave num¬ 
ber k > 0 with E{k) = E^ corresponding to > 0 as 
given by Eq. Apart from this, there is a pair of stag¬ 
gered evanescent modes with wave numbers K 2 = tt ± m 
where k > 0 and E{K 2 ) = E. Generally, the interplay 
of an open and a closed scattering channel is expected 
to lead to resonance effects, as have been observed for 
instance in numerous variations of tight-binding lattices 
with side-coupled defects [53]. In the following, we will 
explicitly work out the general transmission properties of 
the zigzag-defect model and in a second step specialize 
the results to the weak-coupling resonance regime. 

We first note that cos +cos A2 = — ti/(2t2) according 
to Eq. ([^, so the parameters of the travelling and the 
evanescent wave, respectively, at the same energy are re¬ 
lated through 


K, = acosh 



+ cos k 


where k < kc ensures that k is real. To obtain the sta¬ 
tionary scattering solutions, we employ the Ansatz: 


(t)j = , j < 0, 

(pj = j > 0, 


where R, T, Ai , A 2 and 0o are unknown complex num¬ 
bers and the signs have been chosen such that there 
is no asymptotic exponential growth. Inserting this 
into Eqs. § now results in a 5 x 5 inhomogeneous 
linear system of equations which can be solved for 


(R, T, Ai, A 2 , 0o) analytically, yielding 


R 




sink 


+ 


j‘^v{k) 


\sinh/i: VE{'y‘^ — 1) 

T = R+1 

^‘^v{k) 


1 + i 


sin k 


sinh/i: V + ^(7^ — 1) 


n -1 


- 1 ' 


. . . sink ^ 

Ai — A2 — 

Sinn K 

(/)o = (i? + Ai + l)/7. 


(13) 

(14) 

(15) 

(16) 


Here, v{k) denotes the group velocity at k as given by 
Eq. In the limit of ^2 ^ 0 (such that sinh/i: ^ oc), 
we recover the result of [52] for a nearest-neighbor chain 
disturbed by a defect. Taking also 7 = 1, the above ex¬ 
pressions reduce to those for a potential-only (no modula¬ 
tion of the hopping) impurity in a nearest-neighbor chain, 
R = [2itisink/V — 1]~^, T = [1 + W/(2ti sin/c)]~^, 
Ai = A 2 = 0 . 

Using the continuity equation ([^, we find that the con¬ 
tributions of the evanescent modes to the current cancel 
out and we are left with the probability conservation re¬ 
lation |Rp + |Tp = 1, where the probability for reflection 
is explicitly given by 


PR = \R\^ = 


sin k 
sinh R 


^^v{k) 


V + E{j^ - 1) 


-1 


(17) 

We can immediately read off a number of features here. 
Eirst, if U = 0 and 7 = 1 we recover the trivial case of 
the homogeneous lattice with = 0 at all energies (we 
will exclude this case in the following). Eurther, if 7 = 0 
(i.e., the defect site is fully detached from the lattice), the 
transmission and reflection probabilities are independent 
of U, as expected, and interestingly only depend on the 
ratio sink/ sinh/i:. Eor arbitrary U, 7, in the low energy 
limit of /c ^ 0 we find ^ 1, so the transmission prob¬ 
ability Pt = ^ — Pr vanishes, and neither the defect site 
nor the evanescent mode are n otably excited, ^ 0, 
Ai 0, as is seen from Eqs. ( 15|16 ). In contrast, ap¬ 
proaching the upper edge of the energy interval, E ^ Ec^ 
we need to distinguish two cases. Eor ^2 < U/d, Ec co¬ 
incides with the upper band edge and the corresponding 
/c = TT, while R remains finite. Thus, p^^ ^ 1 and the 
transmission probability drops to zero when approach¬ 
ing Ec. In contrast, for ^2 > ti/4 we have sinh/i: ^ 0 
when E approaches Ec from below, while sin kc and v{kc) 
remain finite and thus the transmission probability be¬ 
comes unity. This full transmission effect is independent 
of the detailed features of the defect in that it persists for 
arbitrary U, 7. Comparable transmission resonances at 
energies where a new scattering channel opens have been 
observed in transverse-multimode waveguides [ZIIIUISI], 
see also [551 [sg. The full transmission at Ec is accompa¬ 
nied by Ai ^ — 1, so the evanescent waves (whose decay 
length goes to infinity ns E ^ Ec) are excited at an am¬ 
plitude identical to that of the incoming travelling wave. 
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Furthermore, ^ 0- All these features can be observed 
in Fig. 



Figure 6 . (Color online) Parameters of stationary scattering 
solutions in the subcritical energy range Eq < E < Ec'. Prob¬ 
ability of transmission/reflection, squared amplitude of the 
evanescent mode and occupation of the defect site. Model 
parameters are chosen asti = 1, 7 = 0.5, V = — 2 , and (a) 
t2 = 0.2, (b) t2 = 0.5. 


Apart from the limiting cases at the edges of the energy 
interval, it is seen that the transmission probability also 
becomes unity if V ^( 7 ^ — 1) = 0. Thus, if the energy 
Et = V/{1 — lies in the interval Eq < Et < this 
produces a transmission resonance inside the band. No¬ 
tably, this can only occur for 7 7 ^ 1 . Similarly, = 1 if 
V -i-E(j^ — 1) +7^(2ti + 8^2 cos k) sinh k = 0, leading to a 
reflection resonance inside the band if the energy Ej^ that 
satisfies this equation lies in the interval Eq < Er < Ec. 
In the weak-coupling limit of 7 ^ 0 the latter transmis¬ 
sion and reflection resonances both approach V. Thus, if 
now the defect energy satisfies Eq < V < Ec^ for E ^ V 
the transmission properties of an incoming wave will de¬ 
pend very sensitively on its energy, as is most clearly seen 
in the corresponding wave packet dynamics. So far, we 
have studied stationary scattering solutions. Superim¬ 
posing these, weighted by a distribution that is localized 
in /c-space, immediately gives insight into the dynamics 
of wave packets when scattering from the impurity m- 
In particular, for a wave packet that is well localized in k- 
space (and thus broad in direct space), the reflected and 
transmitted fraction are immediately determined by the 
reflection and transmission probabilities of the station¬ 
ary scattering solution at the central k. We check this 
by initializing wide Gaussian wave packets centered at 
different k near the expected resonance and propagating 
them towards the defect by direct numerical integration 
of the discrete Schrodinger equation 0 . Explicitly, the 
initial condition reads ipj oc exp [—{j — E ikj'\ 

with jo <C 0 the initial central position and the width 
parameter rco ^ 1. The results are shown in Fig. 
demonstrating almost perfect transmission/reflection of 
the wave packet from the defect upon a small variation of 
the incoming energy near the resonance, accompanied by 
a strong transient excitation of the defect site during the 
scattering process. It is in this weak-coupling limit that 
the zigzag-defect model provides a clean example of the 
Fano-Feshbach resonance mechanism. We will analyze 


this in detail in the following. 



0 20 40 60 80 100 


T 

Figure 7. (Color online) Wave packet dynamics at neighbor¬ 
ing transmission and reflection resonances for small 7 , the 
color encodes |'0n|^- (a) k = 0.66, E = —2.076. (b) k = 0.71, 
E = —1.817. Other parameters are ti = t 2 = 1, F = —2, 
7 = 0.2, wq — 100, jo = —300. The insets show snapshots 
of the respective densities at time r = 60, demonstrating the 
resonant excitation of the defect site n = 0 . 


Fano resonance at weak coupling To make the con¬ 
nection to the canonical Fano resonance formalism (see 
e.g. m), it is useful to rewrite the transmission coeffi¬ 
cient as T = F/{F — iG)^ where 

F{E) = f + f;(72 - 1) 

+ 7 ^ ( 8^2 cos/c + 2 ti) sinh/^, (18) 

are real functions, cf. m- Then, the aforementioned 
transmission/reflection resonances Er^Er occuring near 
V at weak coupling correspond to the zeros of G and F, 
respectively: G{Et) = 0 , F{Er) = 0 . To lowest non¬ 
vanishing order in the coupling parameter 7 ^, one can 
expand the zeros as 


Et ^ F + 7V, 
Er^V + 7V 


+ 7^ 8^2 


^cos kv + 




+ cos kv 


- 1 , 

































where 


cmkv = -ji + 


At^J 4t2 2 


In a second step, one can linearize F and G around 
their respective zeros, assuming that these zeros he close 
enough to each other such that there is a common re¬ 
gion where both linearizations apply (which is certainly 
true in the limit of 7 ^ ^ 0 in which both zeros coincide 
and will approximately still hold for small enough 7 ^). 
Then the transmission coefficient in the vicinity of the 
resonances is given by 


T = 


F\Er){E-Er) 


F'{Er){E - Er) - iG'{ET){E - Et) 


within the weak-coupling expansion, where the prime de¬ 
notes a derivative with respect to E. This yields the 
well-known Fano profile 


Pt{E) = \T{E)\^ = - 

upon the identifications 

G'{Et) ^ Er - 
q = „„„ , , E = 


F'iEnY 


1 {E-E + qT/2f 
+7 (r/2)2 + (i; - E)2 


q^Er „ Et — Er 
,T = 2q- 


1 + q 


1 + q^ 


( 20 ) 


( 21 ) 


To evaluate this, we calculate the zeroth order contribu¬ 
tions of the derivatives which we find to be 


F'{Er) ft 

'-1, 

( 22 ) 

G\Et) ft 

tl E Atit2 COS kv 
y 4^2 (1 — cos^ kv) 

(23) 


To lowest order in the coupling, this gives for the Fano 
asymmetry parameter 


Q = 


— At2V + 8^2 — 4^2 — 2Vt2 + t\ 
ti \/t\— At2V + 8^2 + 4^2 + 2 V ^ ^2 — t\ 


1/2 


(24) 

Corresponding approximate values for the resonance po¬ 
sition E and width F can then be obtained by inserting 
this into the expressions of Eq. (21). The ensuing predic¬ 
tion of the Fano resonance profile, Eq (20), is compared 
to the exact transmission line in Eig. for an example 
set of parameters, showing very good agreement. As ex¬ 
pected, in the weak-coupling limit the width F is of order 
7 ^, thus producing the characteristic narrow resonance 
structure. 

Eq. suggests a high degree of tunability of the 
Eano lineshape through the asymptotic hopping parame¬ 
ters ti, ^2 of the lattice and the defect potential V. Eirst, 
we note that in this system q < 0 throughout, thus there 
is no g-reversal even when V crosses zero. The absolute 
value of g, on the other hand, is subject to large vari¬ 
ations as illustrated in Eig. When concentrating on 
small t 2 and |E|, we find the limiting case discussed in 



Figure 8. (Color online) Fano resonance at weak coupling, 
ti — t 2 — 1,7 = 0.1, V = —2. (a) Transmission probabil¬ 
ity pt as a function of the incoming energy in the subcriti- 
cal inter val Eq < E < Ec- (b) Functions E and G defined 
in Eqs. ( 18|19 ) near resonance (solid black lines), together 
with their zeros Er^ Et and linearizations as predicted from 
the weak-coupling expansion, Eqs. ( |22p^ (dashed red lines), 
(c) Zoom to shaded region of subfigure (a) together with the 
predi cted F ano profile within the weak-coupling expansion, 
Eqs. ( 20|21[ ) (dashed red line). 


[l3] in which g <C — 1 and the Eano profile approaches a 
symmetric Breit-Wigner lineshape [Eig. [^a)]. The same 
is found more generally when V ^ Eq and thus the res¬ 
onance occurs at small k. Interestingly, in the opposite 
limit ofV<Ec the limiting value of g crucially depends 
on the ratio of the hopping parameters: If t 2 < ti/4, 
then again g <C — 1 , eventually approaching a symmet¬ 
ric Breit-Wigner lineshape [Eig. |^b)]. But if t 2 > ti/4, 
we find that g ^ 0 instead, asymptoting towards an in¬ 
verted Breit-Wigner shape [Eig. |^c)]. This qualitative 
change in g reflects the qualitative change of the trans¬ 
mission probability at Ec for the two cases (full trans¬ 
mission for t 2 > ti/4, zero transmission for ^2 < ^i/4). 
Eor the former case, the resonance reduces to an approx¬ 
imately symmetric dip on top of a background of trans¬ 
mission unity, while for the latter case it corresponds to 
an approximately symmetric spike of full transmission 
on a zero background. In contrast, a maximally asym¬ 
metric Eano lineshape with g = — 1 is obtained when 
V = tlli2t2)-2t2 [Fig.§d)]. 


VI. SCATTERING AT SUPERCRITICAL 
ENERGIES 


We now consider scattering from the defect for energies 
at which two different propagating waves coexist. As 
discussed before, this can only happen if ^2 > ti/4 and 
Ec < E < Em, which is assumed throughout this section. 
Let ki denote the wave number of the incoming wave and 
k 2 (of the same sign and thus opposite group velocity) the 
second solution with E{ki) = E(k 2 ) = E. The scattering 
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Figure 9. Transmission probability pr as a function of energy 
Eq < E < Ec in the weak-coupling regime, demonstrating 
tunability of the Fano g-parameter via t 2 ^V. Throughout 
ti = l,7 = 0.1,t2 and V as given in the figures. The indicated 
values of q are calculated from Eq. (24) and agree with the 
result from a Fano profile fit to the respective resonance region 
on the percent level. 


Ansatz is now 


(j)j = j < 0 , 

+T2e-*''=^ j > 0, 


where again is an additional unknown. The terms with 
prefactors Ti correspond to refiected/transmitted 
waves in the two channels, respectively, while only 
the exp(i/cij) wave is incoming. Note the different 
choices of signs accordings to the different signs of the 
group velocities. This Ansatz assumes v(ki) > 0 and 
v{k 2 ) < 0 , thus one of the following two cases must 
hold, see Fig. Either ki G {kc^ km) and k 2 = 
acos [—ti/( 2 t 2 ) — cos/ci] G (/c^, 7 r), or ki G {—ir^—km) 
and k 2 = —acos [—ti/{2t2) — cos ki] G {—km^ —kc)- 
Inserting the scattering Ansatz into Eqs. ^ and solving 
the resulting inhomogeneous linear equation now yields 


Ri = 


I _ ^ ■ 

V + E{'y‘^ — 1) sin/c2 

sinki ^ ^ 

' - Ti = i?i + 1, 


-Ri, 


R 2 =T2 = — — — 
sm k 2 

00 = (^1 + i?2 + l)/7- 


Evaluating the continuity equation (|^, we find the fol¬ 
lowing relation reflecting probability conservation 


Pt,i Pr,i 




v{ki) 


\T2? 


){k2) 


v{ki) 


\R2? = 1, (28) 


which can be checked to be satisfied by the above explicit 
expressions. It is worth noting that for the two different 


wave numbers of the same energy the ratio of the group 
velocities is given by the simple formula 


v{k 2 ) _ sinfe 
v{ki) sin/ci’ 


iiE{k^)=E{k2)Ak^\^\k2\. (29) 


Formally, replacing k 2 — Eqs. (13)-(15) can be 

recovered from Eqs. ([^-(27). 


Fig. shows the total transmission probability pt = 
Pt,i + Pt,2 = 1 — Pi? as a function of the incoming en¬ 
ergy E when varying the coupling 7 , while keeping ti, 
t 2 , V fixed. We include here both subcritial energies 
Eo < E < Ec (where there is only one open channel) 
and supercritical energies Ec < E < Em- In the latter 
range, specifying the energy E does not uniquely fix the 
incoming wave number (as per our above discussion of the 
degeneracy of the band structure curve). The two differ¬ 
ent initial quasi-momenta corresponding to the same E 
here can be distinguished by their sign, ki > 0 or ki < 0 . 
First, it can again be seen here that at Ec the transmis¬ 
sion probability is unity. In the weak-coupling regime, we 
clearly observe the Fano resonance structure at an energy 
near f/, as discussed before. Increasing the coupling, the 
transmission and reflection resonances are shifted away 
from V. The reflection resonance Ep moves to larger 
E. Depending on the details of the parameters, it may 
reach Ec at a finite 7 ^ and then cease to exist (this hap¬ 
pens at 7 ^ = 1 — if the latter quantity is positive). 

Interestingly, in the case shown here, with Ec = 0, the 
reflection resonance does not vanish at a finite 7 ^, but in¬ 
stead becomes very narrow and asymptotes to Ec (where, 
as noted before, full transmission is found). Thus, also 
at stronger couplings we identify a region in which the 
transmission probability sensitively depends on the in¬ 
coming energy (similar to the Fano resonance interval at 
weak coupling), namely for energies just below Ec- The 
fate of the transmission resonance Ep = V/{1 — 7 ^) is 
mainly determined by the sign of V. For V < 0, as is the 
case in Fig. it drifts towards the lower band edge Eq 
and disappears there nt = 1 — V/Eq. In cases with 
1 — V/Ec finite and positive, the transmission resonance 
reappears at Ec for larger 7 ^. This does not happen here 
due to F^c = 0 in the example. For 0 < "F < E^ Ep first 
drifts towards Ec with increasing 7 ^ instead. 


In addition to the total transmission probability, 
Fig. also displays the respective contributions of the 
second channel in the reflection and the transmission, re¬ 
spectively, see Eq. (28). These can be simplified to 


PR ,2 
PR 
PT,2 
Pt 



I- Pt a ^ sin ^2 
Pt V 


sin k 



2 


1 


-1 


(30) 

(31) 


Remarkably, the contribution of the second channel to 
the reflected fraction is thus independent of the defect pa¬ 
rameters F, 7 and only depends on the incoming momen¬ 
tum ki and the asymptotic parameters ti, t 2 of the lattice 
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E 



Figure 10. (Color online) (a) Total transmission as a func¬ 
tion of the incoming energy E and the coupling parameter 7 . 
The degenerate regions in the supercritical regime are distin¬ 
guished by the sign of ki, see discussion in the text. Note 
that above the critical energy Ec the energy axis is nonmono¬ 
tonic, increasing towards Em and then decreasing to Ec again. 
(b,c) Branching ratios of reflection/transmission into the sec¬ 
ond channel to the total reflection/transmission. The vertical 
lines Bit E = Em separate the regions of /ci > 0 and ki < 0, 
as in (a). Throughout, = t 2 = 1,1^ = —2. Crosses mark 
the parameter values at which wave packet propagation runs 
are shown in Figs. [7|andpT] 

(which determine the corresponding /C 2 ). In Fig. it is 
clearly seen that Pr^Ipr is independent of 7 (see also 
Fig, below). In contrast, the contribution of the sec¬ 
ond channel to the transmitted wave does depend on the 
overall transmission, and thus on the defect details. 
There are regions of parameter space in which pr^^Ipr or 
Pt, 2 /pt approach 1 , i.e. where the reflected/transmitted 
wave is dominated by the second channel (e.g., most 
prominently for the reflected part at ki > kc). This, 
however, tends to be accompanied by small overall re¬ 
flection/transmission: Only in regions where just a small 
part of the incoming wave is reflected/transmitted will 
the reflected/transmitted part have dominant contribu¬ 
tions from the second channel. For more generic param¬ 
eter values with incoming energies Ec < E < E^^ there 
are nonvanishing contributions of both the first (incom¬ 


ing) and the second channel in both the transmitted and 
the reflected part. This is observed in the wave packet 
propagation runs shown in Fig. [m A wave packet of cen- 




Figure 11. (Color online) Wave packet dynamics in the su¬ 
percritical regime, the color encodes |' 0 n|^- (a) ki = 1 . 2 , 

E = 0.75. (b) ki = —2.7, E = 0.539. Other parameters: 
ti = t 2 = 1, 7 = 1, "F = — 2 , jo = —300, Wo = 25. 

tral momentum ki with E{ki) > Ec is initialized in the 
asymptotic region of the lattice. It propagates towards 
the defect at the group velocity v{ki). After scattering 
from the defect, there are four wave packets, two cor¬ 
responding to the incoming channel and propagating at 
v{ki) (transmitted) and —v{ki) (reflected), respectively. 
In addition, two new wave packets are emitted symmet¬ 
rically from the collision event, equal in shape (accord¬ 
ing to T 2 = R 2 ) and travelling at the velocities ±.v{k 2 ). 
Fig. also displays the final contributions of the four 
individual wave packets to the total norm. These agree 
with the weights in Fourier space, as expected. When the 
initial wave packet is wide enough to be considered es¬ 
sentially localized in momentum space, the weights of the 
individual wave packets after the separation agree well 
with the reflection/transmission probabilities extracted 
from the stationary scattering solutions (within below 
1% for the runs shown in the figure). We find that in 
all cases considered pr ^2 + Pt ,2 <0.5, thus the transfer 
to the second channel does not exceed 50%. A qualita¬ 
tive difference between the regimes with incoming ki > 0 
and < 0 is worth noting. In the former case, we have 
> \'^{k 2 )\^ so the secondary wave packets created 
in the scattering process lag behind those travelling at the 
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velocity of the incoming wave packet [Fig. IT^a)]. In the 
latter case of —it < ki < —km-, we find \v{J^\ > \v{ki)\ 
instead, so the emerging wave packets at momenta ±/c 2 
travel at a velocity faster than the initial one [Fig.pT][b)]. 

Finally, let us have a look at the transmission and reflec- 



Figure 12. (Color online) (a): Total transmission as a func¬ 
tion of the incoming energy E and the defect potential V. 
(b,c) Branching ratios of reflection/transmission into the sec¬ 
ond channel to the total reflection/transmission. Throughout, 
ti = t 2 = 1,7 = 0.8. Crosses mark the parameter values at 
which wave packet propagation runs are shown in Fig. |13| 

tion probabilites as a function of the defect potential V 
at a fixed coupling 7 , as shown in Fig.[^ A remarkable 
symmetry property of the scattering coefficients is found 
if 7 = 1 (no modulation of the hopping amplitudes due 
to the defect), in which case V ^ —V only causes com¬ 
plex conjugation of Ti^2 5 ^1,2 in the E > Ec regime and 
thus the scattering probabilities become independent of 
the sign of the defect potential. The remnants of this 
perfect symmetry at 7 = 1 can be seen in Fig. which 
shows results for 7 = 0.8. The second clear feature to be 
observed in Fig. is the region of large transmittivity 
at energies close to (and especially just below) Ec^ inde¬ 
pendently of the defect potential V, as discussed above 
in Sec. |V| In Fig. we demonstrate that wide wave 
packets localized at momenta ki ^ kc indeed can pass 
the defect with little disturbance for vastly different val¬ 
ues of V. The small deviations from perfect transmission 


we attribute to the inevitable finite width in momentum 
space, together with the fact that for ki > kc the trans¬ 
mission probability rather quickly drops below unity. 




Figure 13. (Color online) Almost full wave packet transmis¬ 
sion when approaching the critical energy, k = /cc — 0.01. The 
color encodes |'0n|^- (a) V = +5. (b) V = —5. Other pa¬ 
rameters are: = t 2 = 1, 7 = 0.8, jo = —300, wq = 100. 


VII. SUMMARY AND CONCLUSIONS 

We have investigated the scattering problem in a dis¬ 
crete Schrodinger lattice with first- and second-neighbor 
hopping and a single-site defect, as can be realized 
with evanescently coupled optical waveguides or ultra¬ 
cold atoms. Although the lattice is multiply connected, 
the continuity equation could be recast into the standard 
one-dimensional form by introducing a suitable general¬ 
ized local current. We explored bound states at the de¬ 
fect, noting in particular the absence of symmetry under 
the staggering transformation. Turning to scattering so¬ 
lutions, an inspection of the band structure of the homo¬ 
geneous system revealed two different energetic regimes, 
separated by a critical energy marked by a transmission 
resonance, independently of the details of the defect. At 
subcritical energies, one closed and one open channel co¬ 
exist and give rise to Fano-Feshbach resonance phenom¬ 
ena at weak coupling, where analytical approximations to 
the Fano lineshape parameters were given that demon¬ 
strated a large degree of tunability of the resonance asym¬ 
metry. At supercritical energies, there are two open chan- 
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nels which are coupled in the scattering process, giving 
rise to peculiar wave packet dynamics, where an incom¬ 
ing wave packet splits into multiple fragments moving at 
different group velocities. The changes of the branching 
ratios of this splitting process when varying the model 
parameters were analyzed. 

In several aspects, our results are reminiscent of those 
found for a cubic two-leg ladder with two transverse 
modes, where longitudinal and transverse motion asymp¬ 
totically decouple m- That system, when perturbed by 
an immersed defect, for instance also exhibits a defect- 
independent transmission resonance at the threshold en¬ 
ergy for the opening of the second channel, and the ana¬ 
log of the wave packet splitting we observe would be a 
partial transfer into the other transverse mode during 
scattering, which redistributes energy from the longitu¬ 
dinal motion (thus providing an intuitive explanation for 
the different group velocities of the outgoing wave pack¬ 
ets). Our model is fundamentally different from this cu¬ 
bic ladder because it does not admit the notion of sep¬ 
arated longitudinal and transverse degrees of freedom, 
yet, to some extent the zigzag lattice may be thought of 
as a continuous deformation of the ladder, cf. [39]. The 
increased connectivity of the one-dimensional chain due 
to the second-neighbor hopping to some extent mimicks 


a second (transverse) dimension. 

We have chosen here the second-neighbor hopping model 
with the single-site defect because it is the most sim¬ 
ple framework to observe the demonstrated effects and 
at the same time admits a transparent analytical treat¬ 
ment. Extensions to other types of defects are straight¬ 
forward, and also a more abstract access in terms of 
Green’s functions has been sketched before gnilSO]. Us¬ 
ing nanofiber-based optical traps of helix-shape for ul¬ 
tracold atoms m, realizations of related discrete models 
with sizable hopping to selected remote neighbors (also 
beyond the second one) may become accessible m, with 
three-dimensional helix arrangements of lattice sites gen¬ 
eralizing the zigzag. In that case, the band structure may 
exhibit multiple extrema within the first Brillouin zone 
and more channels are added to the problem, suggesting 
the existence of multiple scattering resonances and en¬ 
hanced wave packet fragmentation in scattering events. 
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